Calculation of the spectra of the superheavy element Z=112 
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Accurate ab initio calculations of the energy levels of the superheavy elements Z=112 are pre- 
sented. Relativistic Hartree-Fock and configuration interaction methods are combined with the 
many-body perturbation theory to construct the many-electron wave function for valence electrons 
and to include core-valence correlations. Two different approaches in which the element is treated 
as a system with two or twelve external electrons above closed shells are used and compared. Sim- 
ilar calculations for mercury are used to control the accuracy of the calculations. The results are 
compared with other calculations. 



I. INTRODUCTION 

Study of the superheavy elements with nuclear charge 
Z > 100 is an important area of research motivated 
by the search for the island of stability (see, e.g. P, 0, 
yj). Synthesis and investigation of superheavy elements 
are conducted at leading nuclear-physics laboratories in 
Dubna, Berkeley, Darmstadt and others. Elements with 
nuclear charge up to Z=118 have been synthesized (4j. 

The study of the element Uub (Z=112) is an impor- 
tant part of this research. Since it was synthesized in 
Darmstadt in 1996 [![ there were numerous works dis- 
cussing its productioUjiiuclear and chemical properties, 
etc. (see, e.g. Refs. H 0, i, i, M, M, El and refer- 
ences therein). In contrast, only very few works devoted 
to the study of electron structure and optical spectrum 
of the element. Eliav et al [ll calculated ionization po- 
tential of neutral Uub and few low energy levels of Uub+ 
and Uub^+. More detailed study of neutral Uub were 
recently reported in Refs. [13, Eli- Quantum electrody- 
namic corrections (QED) for the Uub element were stud- 
ied in Ref. 

In present paper we try to address the shortage of data 
on the electron structure and energy spectrum of the Uub 
element by calculating its energy levels. Element 112 has 
electron structure similar to those of mercury. Therefore, 
we use the calculations for mercury as a test of the calcu- 
lations and as a guide for their accuracy. Most of lower 
states of both atoms can be considered as states with two 
valence electrons above closed shells. We use the com- 
bined configuration interaction and many- body pertur- 
bation theory method (CH-MBPT) [l^, [lil to perform 
calculations for such states. This method has been suc- 
cessfully used for many different atoms [l^, [13, [2TI.[2^ . [23j 
including the superheavy element with Z=120~[24l|. 

There are also states in mercury and element 112 with 
excitations from the 5d or 6d subshell. They cannot be 
considered as two-electron states and in this case we use 
a version of the configuration interaction (CI) technique 
with has been developed for atoms with open d or f 
shells (25l . [26j . Some states are covered by both methods 
which is another test of the accuracy of the calculations. 
We also compare our results with the calculations of Li 
et al in Ref. 



II. METHOD OF CALCULATIONS AND 
RESULTS FOR MERCURY 



Many states of mercury and element 112 (E112) can be 
considered as having two valence electrons above closed 
shells. The uppermost core subshell is the 5d^° subshell 
for mercury and the 6(i^° subshell for E112. However, 
it is well known that mercury also has states of discrete 
spectrum which have one electron excited from the 5^^*^ 
subshell [13. The lowest such state, the 5d^6s^6p ^Pg 
state is obviously due to the 5d^/2 ~^ excitation. Its 
energy is 68886.60 cm"-'^ which is roughly double of the 
minimal excitation energy (see Table HI)) . It is clear that 
the 6dc,/2 — + 7pi/2 excitations in the El 12 superheavy 
element must be even easier due to larger fine structure. 
Indeed, with fine structure increasing the 6(^5/2 and 7pi/2 
states move towards each other on the energy scale. The 
6(^5/2 state goes up while the 7pi/2 state goes down. This 
means that one should expect of having even more states 
with excitations for the 6d^^ subshell in the discrete spec- 
trum of El 12 than those found in mercury. And these 
states are expected to have lower energies. 

The presence of the states with the d — p excitations 
from the core is a serious complication for the calcu- 
lations. The two-valence-electrons atoms like Ba, Ra, 
E120 [H, [m [13, [Hi can be treated very accurately 
by means of the configuration interaction (CI) tech- 
nique combined with the many-body perturbation theory 
(MBPT) (the CI+MBPT method |l3|). In this method 
the CI technique is used to construct the two-electron 
wave function and to include correlations between two 
valence electrons to all orders via matrix diagonalization. 
The MBPT is used to include the core-valence correla- 
tions. This method docs include the core- valence excita- 
tions but in an approximate way, using the lowest order 
perturbation theory. This might be not very accurate in 
the case when states with the core- valence excitations are 
in the discrete spectrum, like in mercury and E112. 

The aim of present work is to predict the spectrum 
of the El 12 superheavy element. Since it has both 
types of states, with and without excitations from the 
6(i^° subshell, we use two different methods of calcula- 
tions. One is the CI-I-MBPT method for two valence 
electrons [13, [H, [HI [13, [H (method A) and another is 
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the CI method for twelve electrons [H, (method B). 
We demonstrate that unless a two-electron state happens 
to be very close in energy to a state with the excitation 
of the d-electron from the core the CI+MBPT method 
gives remarkably accurate results. The twelve-valence- 
electrons method B is used to find positions of the states 
with the excitations from the core. 



The coefficients Ci as well as two-electron energies are 
found by solving the matrix eigenvalue problem 

{H''^ - E)X ^ 0, (7) 

where Hff = ($,|i7°ff |$,) and X = {ci, C2, . . . , c„}. 

The results of calculations for Hg with method A will 
be discussed in section HTCl 



A. CI for two electrons: method A 

Here we use the CI-I-MBPT method developed in our 
earlier works [13, E, HH, H^] . The calculations are done 
in the V^^^ approximation [l9l | which means that initial 
Hartree-Fock procedure is done for a double ionized ion, 
with two valence electrons removed. 

The effective CI Hamiltonian for a neutral two-electron 
atom is the sum of two single-electron Hamiltonians 
plus an operator representing interaction between valence 
electrons: 



H' 



off 



hi{ri) + hi{r2) + /i2(?'i,f2)- 



(1) 



The single-electron Hamiltonian for a valence electron 
has the form 

hi^ho + ti, (2) 
where Hq is the relativistic Hartree-Fock Hamiltonian: 

/lo = cap + (/? - l)mc2 + V^-^, (3) 

r 

and El is the correlation potential operator which rep- 
resents correlation interaction of a valence electron with 
the core. 

Interaction between valence electrons is the sum of 
Coulomb interaction and correlation correction operator 
S2: 



ri -ra 



(4) 



S2 represents screening of Coulomb interaction between 
valence electrons by core electrons. Wc use the second- 
order MBPT to calculate correlation operators Si and 
The details can be found in our earlier works [H, 



mm- 

Two-electron wave function for the valence electrons 
^I' has a form of expansion over single-determinant wave 
functions 



(5) 



$i are constructed from the single-electron valence basis 
states calculated in the V^^"^ potential 

^i{ri,r2) = -^{ilJa{ri)ipb{r2) ~ ■ipb{ri)ipa{r2))- (6) 



B. CI for twelve electrons: method B 

The method used in this section has been developed in 
our earlier works [2^ [2^ . As for the case of two valence 
electrons the method is based on the CI technique. The 
main differences between method A and method B are in 
the choice of the basis and in the treatment of the core- 
valence correlations (see below). In general, the method 
B is less accurate than method A. However, its strong 
feature is the ability to deal with large number of valence 
electrons. 

The effective Hamiltonian for valence electrons has the 
form 



12 

E 

i=i 
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(8) 



hi{ri) is the one-electron part of the Hamiltonian 

hi=ca-p+(p- l)mc^ h Vcore + SV. (9) 

r 

Here a and p are Dirac matrixes, Vcore is the Hartree- 
Fock potential due to core electrons and SV is the term 
which simulates the effect of the correlations between 
core and valence electrons. It is often called polarization 
potential and has the form 



6V 



2(r4 + a4) ■ 



(10) 



Here ap is polarization of the core and a is a cut-off 
parameter (we use a = ub)- 

The differences between the Hamiltonian ^ in the 
previous section and the Hamiltonian ([5]) are: (a) the 
5d electrons are treated as core electrons in ([1]) and their 
contribution is included into the potential V^~'^ while 
the 5d electrons are treated as valence electrons in (|8]) 
and their contribution is not included into the potential 
Vcore', (b) the E2 operator is not included in ([5]); the Ei 
operator in ([T]) is replaced by a less accurate polarization 
potential SV in 

To construct the many-electron wave function for 
twelve valence electrons we use the Hartree-Fock single- 
electron basis states which are found by the self- 
consistent procedure performed independently for each 
configuration of interest (see Refs. [2^ [2^ for details). 
Table [T] lists all configurations of the valence electrons 
for mercury and El 12 considered in present work. The 
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TABLE I: Even and odd configurations of Hg and E112 and 
effective core polarizability (a.u.) used in the calculations. 



Atom 


oet 


Parity 


Configuration 


Q-p 


ng 


1 


Even 


OO OS 


U.4 




2 


Even 




0.4 




3 


Odd 


5d^°6s6p 


0.386 




4 


Odd 




0.41 


E112 


1 


Even 




0.4 




2 


Even 




0.4 




3 


Odd 


6di°7s7p 


0.386 




4 


Odd 




0.41 



effective core polarizability parameter ap is treated as a 
fitting parameter. Its values for mercury are chosen to 
reproduce the experimental data for energy levels of the 
corresponding configurations. The same values are then 
used for the superheavy element E112. The results for 
mercury will be discussed in next section. 



C. Results for mercury 

Results for mercury are presented in Table ITIl Here ex- 
perimental energies are compared with the energies cal- 
culated within frameworks of methods A and B which 
are described in previous section. Energy levels of mer- 
cury were calculated by many authors before (see, e. 
g. [3, [l^ nil). A review of these calculations goes be- 
yond the scope of present work. In our case mercury 
serves only as a test of the calculations for the super- 
heavy elements Uub. Therefore we included in Tabl e HIl 
the results of calculations of only one other group 
who also calculated the spectrum of Uub (see next sec- 
tion). 

Method A gives very accurate results unless a state of 
interest happens to be very close to another state with 
the same total momentum J and parity and which has a 
hole in the 5d shell. For example, the largest deviation 
of the theory from experiment in method A is for the 
bd}°&s7p state which is close to the ^(P&s^&p ^P^ 
state. These states are strongly mixed, however this mix- 
ing is included in a very approximate way in method A. 
It treats an atom as a two- valence-electron system and 
excitations from the core are included only in the second- 
order of the MBPT in the E operator in the effective CI 
Hamiltonian. Note that this maximum deviation (1102 
cm~^) is only 1.5% of the energy. 

Method B is less accurate, however it gives the posi- 
tions of the energy levels of the states with excitations 
from the 5d subshell which cannot be obtained by method 
A. 

The results of Ref. are closer to our method B 
results. 



III. RESULTS FOR UUB (Z=112) 

The results of calculations for the superheavy element 
Uub (Z=112) a re p resented in Table|TTT]together with the 
results of Ref. [l5|. We also present the g- factors Lande 
in the Table. This includes the calculated g-factors as 
well as g-factors obtained from analytical expressions in 
the LS and jj schemes. The g-factors are useful for the 
identification of the states. The Uub is a superheavy 
element with large relativistic effects. Therefore the jj 
scheme works better for it than the LS one. However, the 
LS scheme is also useful for the comparison with mercury 
for which the LS scheme is commonly used. 

The ^-factors in the LS scheme are given by (non- 
relativistic notations) 



gNFiiJ, L, S) : 



(11) 



L{L + 1) + S{S + 1) 



2J{J+1) 



where L is angular momentum of the atom, S is its spin 
and J is total momentum (J = L -I- S). 

For the case of two electrons the g-factor in the jj 
scheme is given by 

9niJ,Jij2)= (12) 

J{J + i)-.n{j2 + i)+ji{ji + i) 



9NRiji,h, 1/2)- 



+5iV_R(j2,^2,l/2) 



2J(J+ 1) 

■/(■/+1)-Jl(jl + 1)+J2(j2 + 1) 

2J(J+1) 



where ji and j2 are total momentum of each electron 
and J is total momentum of the atom (J = ji -I- j2) and 
gNR is given by pd|) . The formula also works for an 
electron and a hole (e.g, the (6^5/27^1/2)2 state). 

The main difference in the spectra of mercury and Uub 
is due to larger fine structure in the 6d subshell of Uub 
than in the 5d subshell of Hg. This leads to easy excita- 
tion of the 6d^/2 electron and large number of the states 
in the spectrum of Uub which correspond to the Qd^ls^Tp 
configuration. According to calculations in Ref. fisl it 
also lead to the change of the ground state configuration 
of Uub+ as compared to the Hg+ ion. The ground state 
configuration of Uub"*" is shown to be the Qd^ls^ configu- 
ration compared to the hd^^Qs ground state configuration 
of Hg+ . One should also note large negative relativistic 
correction for 7s energy in Uub which makes this state 
to be more tightly bound than the 6s state of Hg. 

The states of the 6d^7s^7p configuration are calculated 
with method B, states of the 6d^°7s8s and 6d^'~'7s8p con- 
figurations are calculated with method A and the states 
of the Qd^^lsTp configuration are calculated with both 
methods. The results of both methods are in good agree- 
ment with each other and in reasonable agreement with 
Ref. [1^. States of the Qd^ls^lp configuration are well 
separated in energy from the states of the same total 
momentum of the Qd^'^lsTp and Qd^^ls^p configurations. 
This means that the mixing between these states is small 
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TABLE 11: Energy levels of Hg (cm^^) 



Present work Other' 



Config. 


Term 


J 


P 6 



















0.000 

















5d^°6s6p 


3po 





37645.080 


37480 


165 


37763 


-699 


38248 






1 


39412.300 


39338 


74 


39442 


-30 


38441 






2 


44042.977 


44287 


-244 


42887 


1156 


49363 


5d^°6s6p 


Ipo 


1 


54068.781 


54263 


-194 


54442 


-373 


57402 


5d"'6s7s 




1 


62350.456 


62181 


169 








5d^°6s7s 







63928.243 


63681 


247 










3po 


2 


68886.60 






70287 


-1400 


70139 


5d^°6s7p 


3po 





69516.66 


69223 


294 












1 


69661.89 


69397 


265 












2 


71207.51 


70106 


1102 








5d^°6s7p 


Ipo 


1 


71295.15 


71213 


82 








5d^°6s6d 




2 


71333.182 


71327 


6 








5d^°6s6d 




1 


71336.164 


71345 


-9 












2 


71396.220 


71383 


13 












3 


71431.311 


71412 


19 












3 


73119.2 






71825 


1294 


71453 


5d^6s'^6p 


Ipo 


1 


78813 






78174 


639 


79357 



"J.G. Li et al, Ref. 0] 
''Experiment, Ref. [27l | 
'^Calculations with method A 

"^A = Eexp - Ea 
'^Calculations with method B 



and should not affect the accuracy of the results. Judg- 
ing by comparison with mercury we expect the results 
for Uub to be accurate within few per cents. 

Accurate calculations for superheavy elements should 
include Breit interaction, quantum electrodynamic 
(QED) corrections and volume isotope shift. However, as 
it has been demonstrated in our previous works 0, [29| 
even for atoms with Z=120 Breit and QED corrections 
are relatively small and extrapolation of the error from 
lighter analogs of the superheavy atoms is likely to pro- 
duce more accurate results than the inclusion of these 
small corrections. The accuracy of pre sent calculations 
is lower than that for Z=120 in Ref. [24. [29| due to com- 
plex electron structure of the Uub element. Therefore, 
these small corrections can be safely neglected on the 
present level of accuracy. This is in agreement with the 
results or Ref. [T(?| in which QED corrections have been 
considered for E112 and found to contribute about 0.5% 



to the ionization potential. 



IV. CONCLUSION 

We have calculated 17 lowest energy levels of the su- 
perheavy element Uub (Z=112). Comparison with simi- 
lar calculations for mercury indicate that the accuracy of 
the calculations is within few per cents. The results can 
be used in the study of the chemical and spectroscopic 
properties of the superheavy element. 
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E 




E 


.y 


E 
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il S^ /o 7Si /o In 
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0.00 





0.00 





0.00 
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1.29 






35785 
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34150 
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1.08 


1.11 
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1.10 


37642 




4 


fedr;/,, 7n3/o)2 


1.25 


1.25 






56131 


1.25 


60366 


6d^"7s7p 


u 


\* '-' \. 1 ' f L / Z J 


0.00 


0.00 


51153 


0.00 


51212 


0.00 
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3po 

^ 1 


1 7.Si /n 7?7i /n 1? 


1.50 


1.33 
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1.41 
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1.33 
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3p2 

^ 2 


1 7Si /O 7j?'^ /O lo 
\ • i/ z ) ' fAj 2 ) 2 
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3po 
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'So 
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0.00 








6d^°7s8p 


3"po 


(7Si/2,8pi/2)o 


0.00 


0.00 


95903 


0.00 










3po 


(7Sl/2,8pi/2)? 


1.50 


1.33 


95084 


1.39 












(7Si/2, 8^3/2)2 


1.50 


1.50 


97342 


1.50 
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